Non-Equilibrium Thermodynamic Description of the Coupling 
between Structural and Entropic Modes in Supercooled Liquids 
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The density response of supercooled glycerol to an impulsive stimulated thermal grating (q=0.63 
(im _1 ) has been studied in the temperature range (T=200-r340 K) where the structure rearrange- 
ment (a-relaxation) and thermal diffusion occur on the same time scale. A strong interaction 
between the two modes occurs giving rise to a dip in the T-dependence of the apparent thermal 
conductivity and a flattening of the apparent a-relaxation time upon cooling. A non-equilibrium 
thermodynamic (NET) model for the long time response has been developed. The model is capable 
to reproduce the experimental data and to explain the observed phenomenology. 

PACS numbers: 05.70.Ln, 78.47,+p, 64.70.Pf. 



Whenever an inhomogeneous temperature or pressure 
field exists inside a substance, heat and momentum will 
flow giving rise to processes, like thermal diffusion and 
sound propagation, which drive the system toward homo- 
geneity. In a normal liquid, at low enough wave-vector 
q (i.e. at typical values of light scattering experiments), 
the time-scales of the two processes are well separated, so 
that sound propagation is adiabatic and thermal diffusion 
is isobaric. Every other microscopic process evolves on 
such a fast time-scale that it enters the dynamic equa- 
tions simply by determining the actual values of ther- 
modynamic derivatives and transport coefficients. The 
situation changes when a liquid is supercooled below its 
melting temperature and the structural relaxation time, 
r a , rapidly grows up upon cooling. When r Q becomes 
of the order of magnitude of the sound wave period, we 
observe phenomena as the sound velocity dispersion and 
sound absorption. These phenomena have been widely 
investigated by ultrasonic and Brillouin spectroscopies 
and commonly described in terms of a relaxing bulk mod- 
ulus or viscosity Ip], [|. Upon further cooling, r a reaches 
the time scale of the thermal diffusion giving rise to a 
complex frequency dependent heat diffusion that is ob- 
served experimentally by specific heat spectroscopy || 
and forced Rayleigh scattering H Though both vis- 
cosity and specific heat relaxations are manifestation of 
the same microscopic process, there is no commonly ac- 
cepted formulation of the dynamic equations in the re- 
gion where structural relaxation and thermal diffusion 
occur on the same time-scale. The main difficulty arises 
when more than one single thermodynamic derivative 
has to be generalized to have a frequency dependence. 
To this respect NET provides a more fundamental ap- 
proach, compared to generalized hydrodynamics, since, 
once the equation of state is properly written in a suitable 
extended parameter space, the frequency dependence of 
generalized thermodynamic derivatives comes out natu- 
rally. Moreover a thermodynamic approach could, hope- 
fully, be related to that thermodynamic picture of su- 
percooled and glassy state which has recently been the 



subject of great theoretical and computational efforts ||. 
Unfortunately whether local thermodynamic equilibrium 
is still meaningfull in the supercooled regime and what 
are the crucial steps in extending the thermodynamic pa- 
rameters space are still open questions. 

In this work we present an Impulsive Stimulated Ther- 
mal Scattering (ISTS) (7| study of liquid and supercooled 
glycerol in a temperature range that covers the region 
where the characteristic times of structural and entropic 
modes become similar. A strong interaction between the 
two modes occurs and we have observed an apparent slow 
down of the steep increase of structural relaxation upon 
cooling, together with a marked non exponential decay at 
long times where the ISTS signal is usually governed by 
thermal diffusion. The NET model we propose is based 
on local thermodynamic equilibrium in an extended pa- 
rameter space j^] . Using literature data, we reproduced 
accurately the experimental ISTS responses at long times 
and explained unexpected features such as the dip in the 
T-dependence of the thermal conductivity of various su- 
percooled liquids |I| |[ ||. 

Glycerol (99.5+%, <0.1% water, Fluka, glass transi- 
tion temperature 185 K) was transferred under nitrogen 
into a Teflon-coated cell with movable windows ITS] . The 
cell was mounted on the cold finger of a cryostat and out- 
fitted with resistive heaters. A platinum resistance ther- 
mometer was immersed in the sample and the tempera- 
ture was kept stable within ±0.1 K. In the present ISTS 
experiment two infrared (A = 1064 nm) short (~ 100 
ps) laser pulses cross each other in the sample volume 
at an angle of ^6° (q = 0.63^m _1 ) and their interfer- 
ence produces an impulsive, spatially modulated, heat- 
ing. The amplitude of the resulting density grating is 
probed by a third CW laser beam (A = 532 nm) im- 
pinging on the induced grating at the Bragg angle. The 
intensity of the diffracted beam as a function of time is 
stored by a digital oscilloscope and averaged over many 
(~ 5000) pulses. The measured signal can be fully as- 
cribed to density dynamics excited through heating since 
electrostrictive and birefringence effects are negligible in 
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FIG. 1: ISTS data of supercooled glycerol at selected tem- 
peratures are reported in logarithmic scale and normalized to 
unity at about t = O.l/is. The inserts show a blow up of the 
short time region, where the structural relaxation gives rise 
to a rising signal. 



glycerol jLlJ [l2| [h3|. Further experimental details are 
reported in Rcf. ||]. ISTS data were collected in the 
temperature range T=200-^340 K. Signals as long as 1.6 
ms were recorded with a time resolution of 1 ns. In the 
present paper we focus on the long time part of ISTS 
signal (t > 0.1/is) where the acoustic transient is over 
and structural and entropy modes evolve isobarically. 
Selected ISTS data are reported in Fig. [j] showing two 
time regions: "short times" (t < lOfis) in the inserts and 
longer times in tha main panels. It can be noted that at 
short times the amplitude of the density grating (ISTS 
signal) increases with a stretched exponential law, due 
to structural relaxation. On lowering the temperature, 
the characteristic time of this rising component stops to 
grow while its intensity vanishes. At longer times the 
heat diffuses and the density grating decays to zero. This 
long time decay, which in absence of coupling is exponen- 
tial (due to the diffusive character of heat equation) , 
now splits into two components. The faster component 
is nearly exponential and its time constant goes through 
a maximum and gives rise to a dip in the apparent ther- 
mal diffusivity, as already reported for OTP in ||. The 
slower component shows a strong non-exponentiality. It 
flattens out and its intensity disappears as the tempera- 
ture is lowered. To be more quantitative, we fitted the 
data with two stretched exponentials for the " structural" 
rising component and the long time tail, and a simple ex- 
ponential for the intermediate time component (apparent 
thermal diffusion). The decay rates [|l4| of the short and 
intermediate components are reported in Fig. ^. At high 
temperatures, the structural relaxation decay rate (o) is 
described by a Vogel-Tamman-Fulcher (VTF) law (full 
line). The parameters 5=2260 K, T VTF =131 K of VTF 
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FIG. 2: Temperature dependence of charachteristic rates for 
the short rising component -apparent structural relaxation- 
(o) and intermediate exponential decay -apparent ther- 
mal diffusion- (•). Solid line is the VTF law 7 VTF = 
To 1 exp[~2260K/(T - 131K)] 

law, j VTF = Tq 1 exp[— B / (T — Tvtf)} are taken from di- 
electric spectroscopy [[l5| while to — 1-4- 10~ 15 s is scaled 
to fit the data. At temperatures lower than ~240 K the 
short time rising component ceases to represent the struc- 
tural relaxation and flattens around a value of 4 /is. At 
about the same temperature, the apparent thermal decay 
rate deviates from the expected smooth behavior (see the 
insert of Fig. ||) and exhibits a dip at T «230 K. The over- 
all scenario depicted in Fig. [2] suggests the existence of 
an interaction between the structural and thermal relax- 
ation dynamics that leads to a complex relaxation time 
pattern. In the following we will introduce a NET model 
that will allow us to compute the temperature evolution 
of this relaxation times pattern and find a very good 
agreement with the experimental observations. NET || 
provides a very powerful framework to study irreversible 
processes such as heat conduction, diffusion and viscous 
flow, from an unified point of view. However, in this 
formalism, it is not straightforward to consider the non- 
exponentiality observed in the structural relaxation dy- 
namics. On the contrary, the presence of a large number 
of internal relaxing variables can be easily taken into ac- 
count provided that local thermodynamic equilibrium is 
valid in the extended parameter space. Therefore, follow- 
ing Allain et al.[[l6|, we choose to represent the observed 
non-exponentiality as the result of the superposition of 
N linearly relaxing variables jl7|] . In this hypothesis the 
Gibbs free energy law per unit mass reads: 

dg = vdp - sdT - E^^df (1) 

where p is the pressure, s the entropy, A 1 the affin- 
ity of the i th relaxation process and £ 8 the correspon- 
dent progress variable (or order parameter). As noted 
above, we are interested in the time region where the 
pressure became and stays uniform. In this time region, 
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the linearized NET equations M written in terms of the 
q-components of the thermodynamic variables (e.g. T 
stands for T(t) = J exp(iqr)T(r,t)), are: 
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T oPo (ds/dt) 
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2 XT 
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(2) 



where To (po) is the average temperature (density), A the 
thermal conductivity and f3 l phenomenologic constants. 
The first of (||) comes from the linearized combined mass 
momentum conservation laws |l9| . The second and third 
equations of (||) represent the energy conservation law 
and the phenomenological relations for the relaxation 
processes respectively. In order to close the above set 
of equations, we use the local thermal equilibrium in the 
extended parameter space: 
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For simplicity, we assume that the thermodynamic affin- 
ity A 1 does not depend on £ J for j ^ i. Differentiating 
the above equations and substituting in (H) we obtain: 



P = 
dT/dt = 
dC/dt = 



-p a~T + (p 2 oC-/T )Y / l 1 ^(C P /CT)m 
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where we have introduced the following symbols: 
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The ISTS density response is obtained by solving (|5|jg) 
for the initial condition T(Q) ^ 0, £ l (0)=0 and then sub- 
stituting in (^). To reduce the number of parameters, we 
use the simplifying assumption that £, p /^t IS independent 
of i In that case, p{t) can be written as 



p(t)/p(0) = -T(t)/T(0)- A,i 



Eti^C(t) (8) 



a°° Ac, 



where Aa and Ac p are the jump from the relaxed to 
the unrelaxed (with respect to the values of the cor- 
responding thermodynamic derivatives. It can be easily 
shown that in the limit F l R <§; - i.e. at low T or high q 
values- this model predicts a density-density correlation 
function which at long times decays as: 



■ 9 W«£tiA l exp[-ry] 



(9) 



We know from photocorrelation experiments, Mode Cou- 
pling Theory and molecular dynamics simulations, that 
the above correlator is very well described by a stretched 
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FIG. 3: ISTS data (o) and the predicted signal (solid line) 
for four different temperatures. The agreement is very good 
through the whole structural-entropic coupling region 



exponential exp[— {t/T a )^\ where (3 is slightly changing 
with temperature and r a obeys a VTF law. Relying on 
this considerations, we arbitrarily choose a distribution 
of 7V=I50 logarithmically spaced rates r^rjlO**, 
Xi= -I+i/28, with weights A* such that the sum (M) re- 
constructs a stretched exponential with (3=0.65 [pT PI- 
This determines the weights apart from a constant fac- 
tor which in turn can be easily fixed by the value of 
Ac p /c~ = A 1 form specific heat spectroscopy data[[|. 
The temperature simply changes the value of T Q which 
is assumed to obey the already quoted VTF law. We 
can assume (T)=(c° p /c^)T° H (T) where T° H (T) is the 
extrapolation to the whole temperature range of ISTS 
thermal decay rates at high temperature. Finally, from 
p(T) data across the glass transition temperature j2lj one 
finds Aa ~ 3.2a°° . We are now left with no more free pa- 
rameters: for each temperature we can compute the ISTS 
signal and compare it to the experimental data. As exam- 
ples, the results of this comparison are shown in Fig. || for 
four different temperatures, showing an excellent good 
agreement between the data (o) and the model (solid line) 
in the whole examined temperature range. A comple- 
mentary and more insightful way of representing complex 
time responses consists in performing an inverse Laplace 
transform analysis: I(t) = G(log7) exp[— jt] dlogj. 
In other words, one can think of the ISTS signal as a 
superposition of exponentials and ask how the weight 
function G evolves across the coupling region. In this 
representation it is easier to visualize the analyzed phe- 
nomenology as the interaction between an exponential 
thermal process and a broad distribution of relaxing vari- 
ables. As a consequence, the relation between the flatten- 
ing in the "structural" rate, the dip in the thermal rate 
and the appearence of a long time tail becomes evident. 
Such an analysis has been carried out on the simulated 
signals and is summarized in Figs. The top figure in the 
right panel represents G in the high temperature region 
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FIG. 4: Computed evolution of the distribution of rates G as 
a function of the logarithmic ratio r = log 10 (r a /Tf^) between 
structural relaxation characteristic rate and the infinite fre- 
quency thermal decay rate. Circles in the left panel represent 
the experimental values 

where no coupling is present. Structural relaxation man- 
ifests itself as a rising (negative weight) stretched (broad 
distribution) exponential, while at longer times (smaller 
rates) thermal diffusion contributes to the signal with an 
exponentially (narrow) decaying (positive weight) com- 
ponent. As the temperature is lowered, the broad struc- 
tural mode moves to shorter rates until the tail of its rate 
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In conclusion, a NET model, based on the assump- 
tion of local thermodynamic equilibrium in an extended 
parameter space, accounts for the rich phenomenology 
observed in the ISTS experiments. In particular, using 
new ISTS data on supercooled glycerol in the tempera- 
ture region where the structural relaxation and the ther- 
mal diffusion process take place on the same time scale, 
we have demonstrated that the model is able to repro- 
duce the experimental data using literature data from 
other experiments. Further investigations on the possi- 
bility of assuming local thermodynamic equilibrium in 
supercooled liquids are crucial in the development of a 
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